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JS  SUMMARY 

A  flight  investigation  aimed  mainly  at  determining  the  pitch  damping  of  a  60  degree 
Delta  wing  aircraft  between  M  =  0-7  and  M  =  1-2  has  been  undertaken.  Flight  tests 
were  conducted  by  ARDU  and  the  measurements  (Attained  have  been  successfully  analysed 
at  ARL  using  a  modified  Newton-Raphscn  parameter  estimation  procedure.  The  results 
for  pitch  damping  as  well  as  other  longitudinal  aerodynamic  derivatives  have  been  compared 
with  wind  tunnel  and/or  theoretical  estimates.  Further  information  has  also  been  gained 
relating  to  the  necessary  conditions  for  successful  extraction  of  parameters  from  flight  data. 
The  parameter  estimation  procedure  has  proved  to  be  an  effective  analysis  method  even 
when  only  a  modest  amount  of  flight  data,  obtained  using  a  relatively  simple  instrumentation 
package,  has  been  available. 
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1.  INTRODUCTION 


In  an  earlier  note1  dealing  with  the  analysis  of  longitudinal  flight  test  data  for  a  delta  wing 
aircraft  at  Mach  numbers  of  0*71  and  0-96,  a  substantial  discrepancy  was  identified  in  pitch 
damping  at  M  =  0-71  between  the  value  extracted  from  flight  and  that  previously  accepted. 
As  the  results  at  M  =  0-96  were  in  reasonable  agreement,  the  variation  of  pitch  damping  with 
Mach  number  over  the  transonic  speed  range  was  in  question.  The  pitch  damping  derivative 
contributes  both  to  the  transient  response  of  an  aircraft  to  control  input  and  also  to  the  steady 
state  manoeuvring  performance  such  as  a  turn  or  pull-up,  and  is  consequently  an  important 
part  of  any  mathematical  representation  of  an  aircraft.  Since  it  is  difficult  to  measure  in  the  wind 
tunnel  or  to  calculate  theoretically,  especially  at  transonic  speeds,  estimates  for  use  in  a 
mathematical  model  need  to  be  validated  against  flight  data.  In  order  to  do  this,  a  flight  test 
program  was  initiated  to  measure  the  pitch  damping  over  the  Mach  number  range  M  =  0-7 
to  1-2.  The  flight  tests,  including  instrumentation  and  data  recording,  were  conducted  by  the 
RAAF  Aircraft  Research  and  Development  Unit  (ARDU)  at  Salisbury,  and  the  data  obtained 
were  analysed  at  ARL.  The  results  of  these  tests,  which  produced  data  for  other  longitudinal 
derivatives  in  addition  to  pitch  damping,  are  reported  here. 

The  flight  data  were  analysed  using  the  modified  Newton-Raphson  parameter  estimation 
procedure  first  documented  in  Reference  2.  This  is  a  response  error  method  which  compares 
model  output  with  measured  responses  and  systematically  adjusts  the  model  parameters  in 
order  to  obtain  a  best  fit.  The  parameters  are  assigned  weighted  a  priori  values,  based,  for  example, 
on  wind  tunnel  data,  so  that  the  estimated  values  of  the  parameters  will  tend  towards  their 
a  priori  values  unless  the  flight  data  contains  information  requiring  a  different  value  for  best  fit. 
The  method  is  now  widely  used  in  the  analysis  of  aircraft  flight  dynamic  tests  and  has  been 
applied  at  ARL  to  both  longitudinal1  and  lateral3  flight  data  including  the  analysis  of  fin  loads.3 
Reference  4  has  shown  that  the  method  can  also  be  used  to  obtain  lift  and  drag  characteristics 
from  dynamic  flight  manoeuvres.  A  more  detailed  summary  of  the  procedure  as  applied  in  the 
present  investigation  can  be  found  in  Reference  1 . 

Section  2  of  this  Note  describes  details  of  the  mathematical  model  used  in  the  analysis. 
The  model  includes  provision  for  non-linear  variation  of  pitching  moment  and  lift  curves  with 
incidence.  Section  3  deals  with  the  test  data  details,  including  instrumentation  and  recording, 
and  describes  the  pre-processing  undertaken  in  preparation  for  the  identification  program.  A 
priori  values  and  weights  used  in  the  matching  process  are  also  summarised.  The  results  are 
presented  and  discussed  in  Section  4  and  comparisons  made  with  wind  tunnel  and  theoretical 
estimates  where  possible. 

2.  MATHEMATICAL  MODEL 

In  this  section  a  mathematical  model  is  developed  which  describes  the  aircraft  longitudinal 
short  period  response  to  an  elevator  input.  A  basic  linear  model  is  extended  to  include  possible 
non-linearities  with  incidence  in  the  pitching  moment  and  lift  curves.  A  further  extension  is 
required  to  deal  with  a  time  lag  in  the  pitch  rate  measurements. 

2.1  Basic  Model 

It  is  assumed  that  the  aircraft  is  trimmed  in  steady  level  flight  with  speed  K«,  incidence  a* 
and  attitude  B,.  Small  disturbance  motions  about  this  initial  state  are  produced  by  elevator 
inputs.  The  linearised  small  disturbance  equations,  in  body  axes,  describing  the  short  period 
response  have  been  derived  in  detail  in  Reference  I.  The  model  state  equations  take  the  basic 
form: 


x  =  i4x+Ba 


0) 


1 


and  the  calculated  response  is: 


y  =  Fx-f  Gu+b 


(2) 


where 

x  is  the  vector  of  state  variables, 
u  is  the  vector  of  control  variables, 
y  is  the  response  vector, 

b  is  a  measurement  bias  vector  taken  as  zero  initially, 

A,  B,  F  and  G  are  matrices  of  parameters  defining  the  model. 

The  parameter  vector,  c,  contains  some  or  all  of  the  elements  of  A,  B,  F  and  G,  the  elements 
of  b  and  the  initial  conditions. 


2.2  Non-linear  Aerodynamics 

The  incidence  ranges  covered  at  each  Mach  number  in  the  present  series  of  tests  are  sum¬ 
marised  in  Figure  1.  Wind  tunnel  tests5  indicate  significant  non-linear  variation  of  pitching 
moment  with  angle  of  incidence  for  the  delta  wing  configuration  under  consideration,  particu¬ 
larly  at  Mach  numbers  near  1.  Following  the  approach  successfully  applied  in  Reference  1, 
the  moment  increment  due  to  incidence  is  modelled  by  a  cubic  equation : 

AM  =  Ma  a  +  M.i  ot2  +  M„«  a?  (3) 

To  a  lesser  extent,  the  Z-force  variation  with  incidence  is  also  non-linear  and  a  cubic 
variation  is  postulated  as  follows: 

A  Z  =  Z«  «  +  Z«.  «2  +  Z„»  a3  (4) 

In  order  to  retain  a  linear  mathematical  model  the  a2  and  a3  terms  are  treated  as  known 
inputs  and  become  part  of  the  input  vector.  Their  values  are  either  calculated  directly  from 
measured  values  of  incidence,  am  or  if  measured  values  of  incidence  are  not  available,  initial 
linear  model  calculations  may  be  used  to  provide  suitable  values. 

In  addition  to  pitching  moment  and  Z-force  non-linearities  it  may  be  speculated,  for  example, 
that  pitch  damping  may  also  be  a  function  of  incidence  requiring  an  equation  similar  to 
Equations  (3)  and  (4)  above  to  model  it.  Each  such  extension  further  complicates  the  basic  model 
and  adds  extra  parameters  to  be  identified.  It  is  clearly  desirable,  from  a  computational  point  of 
view,  to  adopt  as  simple  a  model  as  possible  so  long  as  it  adequately  describes  the  physical 
system.  The  question  of  model  structure  determination  is  an  important  part  of  systems  identi¬ 
fication  addressed  by  a  number  of  authors.*-7'2  Methods  have  been  developed  to  select  a  best 
model  from  a  large  number  of  candidate  models.  In  the  present  case  the  number  of  options  is 
relatively  small  and  can  be  examined  one  by  one.  If  a  significant  improvement  is  obtained  in 
matching  the  test  data  then  the  extension  to  the  basic  model  should  be  retained.  Thus  it  was 
found  that  the  modification  described  by  Equation  (3)  led  to  significant  improvements  in  the 
subsonic  and  lower  transonic  Mach  number  range  and  was  thus  included  in  the  model.  The 
modification  described  by  Equation  (4)  was  much  less  significant  and  it  was  subsequently 
decided  to  retain  only  the  first  two  terms  of  Equation  (4),  On  the  other  hand,  modelling  of 
pitch  damping  as  a  non-linear  function  of  incidence  produced  little  improvement  and  was  thus 
excluded  in  the  interests  of  model  simplicity. 


13  Pitch-rate  Measurement  Lag 

Initial  attempts  to  match  normal  accelerometer  and  pitch  rate  gyro  measurements  focused 
attention  on  an  apparent  lag  of  approximately  0- 1  second  in  the  pitch  rate  measurements  when 
compared  to  initial  mathematical  model  results.  Further  investigation  pointed  to  the  pitch  rate 
gyro  demodulator  unit  as  being  the  source  of  this  lag.  To  account  for  this  the  model  was  augmented 
by  a  first  order  lag  equation  as  follows: 


2 


dqn  ,  4b  G 

dt  +  r  =rq 


(5) 


where  r  is  the  lag  time  constant  and  the  gain  term  G  is  included  to  account  for  possible  calibration 
and/or  gyro  misalignment  errors  in  the  pitch  rate  measurements,  qm.  Both  r  and  G  can  be  included 
in  the  parameter  vector  to  be  identified  if  required.  Greatly  improved  results  were  obtained  when 
the  lag  was  accounted  for  by  Equation  (S)  with  r  initially  set  to  0- 1  and  G  to  I  -0. 

The  complete  model  used  in  the  present  investigation  can  now  be  summarised. 


The  state  equations  are: 
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The  response  equations  are: 
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The  parameter  vector  is: 

c  =  [ZJmVt,  Z'lmV,,  Z,JmVt,  MJIy,  (M'+MJHy,  MtJly, 

G/r,  -  l/r,  Zo,  Mo,  Z'*lmV"  M.>/Iy,  A/..//y]T  (8) 


Z0  and  Mo  have  been  introduced  to  account  for  possible  state  bias.  The  Z«  term  in 
Equation  (7)  was  found  to  make  very  little  contribution  to  the  results  and  was  subsequently 
set  to  zero  in  Equation  (7)  and  omitted  from  the  parameter  vector.  Equation  (8).  The  measure¬ 
ment  bias  vector,  b,  has  been  assumed  to  be  initially  zero  in  Equation  (7)  but  the  biases  in  the 
measurements  of  qn  and  gnt/Vt  were  generally  included  among  the  parameters  to  be  identified 
in  Equation  (8). 


3.  FLIGHT  DATA 

The  form  of  the  data  and  the  test  conditions  are  first  described.  This  is  followed  by  a  discus¬ 
sion  of  the  preliminary  processing  and  corrections  applied  to  the  data.  Finally,  the  parameter 
a  priori  values  and  weights  are  summarised. 


3.1  Test  CeadMias 

Tests  were  performed  at  the  nominal  test  points  summarised  in  Table  I. 
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TABLE  1 

Nominal  Test  Points 


Height  (ft.) 

Mach  Number 

30,000 

10,000 

0-7,  0-8,  0-9,0-92,  0-94,  0-96,  0-98,  1  05,  1  1,  1-2 
0-7,  0-8,  0-9,  0-92,  0-94,  0-96,0-98,  105,  1  1,  1-2 

The  test  flights,  conducted  by  ARDU  pilots,  took  place  in  clear  atmospheric  conditions 
with  a  minimum  of  turbulence.  Two  flights  were  made  at  30,000  ft  and  one  at  10,000  ft.  In 
each  case  the  set  of  Mach  numbers  was  flown  in  the  sequence  shown  in  Table  1 .  The  pilot  was 
asked  to  trim  the  aircraft  at  each  point  and  then,  with  pitch  dampers  off,  apply  a  rapid  stick 
doublet  input  sufficient  to  produce  a  transient  normal  acceleration  increment  of  at  least  ±  1  g. 
In  some  cases  the  pilot  was  able  to  apply  a  second  doublet  input  after  the  response  to  the  first 
had  died  down,  at  the  same  time  checking  that  the  Mach  number  remained  as  constant  as 
possible. 

At  each  test  point  the  pilot  recorded  the  values  of  Indicated  Airspeed  (IAS),  Indicated  Mach 
Number  (IMN),  height  and  fuel  gone  prior  to  initiating  the  trace  recorder.  The  quantities  recorded 
on  the  paper  trace  recorder  included  IAS,  height,  normal  acceleration,  pitch  rate,  port  and 
starboard  elevator  positions  and  pitch  damper  position.  Angle  of  incidence  measurements  were 
not  available.  The  fast  (about  4  inches/second)  recording  speed  was  used  in  most  cases  although 
some  runs  were  recorded  at  a  slower  speed  (about  1  i  inches/second).  All  instrumentation  and 
data  recording  was  undertaken  by  ARDU  personnel.  The  paper  trace  records  were  then  read 
by  Mathematical  Services  Group  at  Defence  Research  Centre,  Salisbury,  and  the  final  results, 
in  Engineering  units,  provided  to  ARL  in  digitised  form  on  magnetic  tape.  It  was  possible  for 
DRCS  to  read  the  traces  at  the  required  sampling  rate  of  40  per  second  even  at  the  lower  trace 
recording  speed. 

The  subsequent  identification  was  based  on  matching  the  time  histories  of  normal  acceleration 
and  pitch  rate  responses  to  the  given  elevator  control  input. 


3.2  Pre-processing 

True  airspeed  (KT)  was  obtained  using  standard  charts  from  IAS  corrected  for  position  error. 
The  corresponding  true  Mach  number,  Mi,  agreed  with  corrected  IMN  readings  to  within 
0  01  for  most  cases.  The  variation  of  aircraft  centre  of  gravity,  mass  and  moment  of  inertia  in 
pitch  were  calculated  from  the  fuel  gone  readings  using  the  mathematical  model  described  in 
Reference  9.  Using  the  same  model,  values  of  aircraft  trim  attitude,  Be,  and  incidence,  *»,  for 
each  test  point  were  established  for  subsequent  use.  In  addition  the  following  pre-conditioning 
was  applied  to  the  measured  data: 

(a)  Normal  accelerometer  readings  were  corrected  for  location  error  according  to  the 
equation : 

m  =  ntm  +  (■*»<?  +  q2)!g  (9) 

where  n,m  is  the  measured  acceleration,  nt  is  the  value  at  the  centre  of  gravity  and  x», 
r»  are  the  co-ordinates  of  the  accelerometer  relative  to  the  aircraft  body  axes.  Measure¬ 
ments  of  q  were  not  available  so  that  the  correction  x%4  was  only  applied  in  the  second 
stage  of  a  two-part  calculation  procedure  as  described  in  Section  4. 1 . 

( b )  Non-zero  trim  values  of  normal  acceleration  and  port  and  starboard  elevator  positions 
were  subtracted  from  the  respective  measured  records  so  that  the  incremental  values 
could  be  compared  directly  with  the  small  disturbance  model  calculations.  In  addition 
the  resulting  port  and  starboard  elevator  settings  were  averaged  to  provide  the  required 
control  input. 
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3.3  A  Priori  Values  and  Weights 

The  a  priori  values  assigned  to  the  aerodynamic  derivatives  in  the  mathematical  model 
were  mainly  based  on  wind  tunnel  measurements  (e.g.  Ref.  S)  or  else  estimated  from  data  sheet 
sources.  The  non-dimensional  values  of  the  derivatives  used  are  summarised  in  Table  2. 


TABLE  2 

A  Priori  Aerodynamic  Derivative  Values  and  Weights 


M 

-c*« 

per  degree 

-c*„ 

per  degree 

—Cma 
per  degree 

-Cm,, 
per  degree 

— (CmB±Cmi) 
per  radian 

0-7 

0  043 

00176 

0  0024 

0  0079 

2-80 

±0  009 

±0  0035 

±0  0005 

±0-0016 

±0-56 

0  8 

0  045 

00187 

0  0030 

0-0081 

2-82 

±0  009 

±0  0037 

±0  0006 

±0-0016 

±0  56 

0-9 

0  049 

00199 

0  0040 

0-0096 

3-27 

±0010 

±0  0040 

±0  0008 

±0  0019 

±0-65 

0-92 

0  050 

0  0200 

0  0043 

0-0098 

2-86 

±0010 

±0  0040 

±00008 

±0-0020 

±0-57 

0  94 

0  052 

0  0202 

0  0052 

0-0101 

2-45 

±0010 

±0  0040 

±0  0010 

±0-0020 

±0-49 

0-96 

0  054 

0  0200 

0  0061 

0-0103 

2-00 

±0  011 

±0  0040 

±0  0012 

±0-0021 

±0-40 

0-98 

0  056 

00193 

0  0071 

0  0103 

1-53 

±0011 

±0  0039 

±0  0014 

±0-0021 

±0  31 

1-05 

0  055 

00165 

0-0105 

0-0099 

1  -98 

±0  011 

±0  0033 

±0  0021 

±0-0020 

±0-40 

MO 

0  055 

00145 

0-0104 

0-0090 

2-02 

±0011 

±0  0029 

±0-0021 

±0-0018 

±0-40 

1  -20 

0  052 

00115 

0-0093 

0  0071 

I  -88 

±0011 

±0  0023 

±0-0019 

±0-0014 

±0-38 

The  static  pitching  moment  derivatives  in  Table  2  are  based  on  a  centre  of  gravity  reference 
position  at  52%  of  centre  line  chord.  For  each  test  point  these  values  were  adjusted  to  the  correct 
centre  of  gravity  position  which  varied  gradually  as  fuel  was  used  from  about  51  -3%  of  centre 
line  chord  at  the  lowest  Mach  number  to  about  49-2%  at  the  highest. 

Of  the  other  parameters  listed  in  Equation  (8)  the  a  priori  values  of  (J/r  and  —  l/r  were  set 
to  10  and  - 10  respectively  while  all  the  other  parameters  being  estimated  were  set  to  zero 
initially. 

The  a  priori  weighting  assigned  to  each  parameter  was  set  by  assuming  a  standard  deviation, 
<t,  equal  to  20%  of  the  parameter  value.  These  are  also  indicated  in  Table  2.  The  relevant  weighting 
was  then  calculated  as  1  /**.  For  those  parameters  initially  set  to  zero  a  very  large  standard 
deviation  was  assumed  leading  to  an  effectively  zero  a  priori  weighting.  For  those  parameters 
which  have  been  assigned  a  non-zero  a  priori  weighting,  a  comparison  can  be  made  with  the 
"Cramer-Rao"  bound  which  is  calculated  by  the  modified  Newton-Raphson  program.*  A 
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"Cramer-Rao”  bound  significantly  less  than  the  a  priori  standard  deviation  leads  to  increased 
confidence  in  the  extracted  value  of  the  corresponding  parameter. 

Finally,  it  is  necessary  to  assign  a  weighting  to  each  of  the  time  histories  being  matched. 
This  reflects  the  measurement  noise  in  the  respective  channels.  In  the  present  case  no  differen¬ 
tiation  was  made  between  measurements  of  the  pitch  rate,  q,  and  normal  acceleration,  gnz/Ve , 
channels  and,  following  Reference  I,  equal  weight  based  on  a  standard  deviation  of  5°„  of  the 
maximum  value  of  each  measured  record,  was  assigned  to  both.  For  records  not  being  matched 
(e.g.  a,  6)  effectively  zero  weighting  was  assigned. 

4.  RESULTS 

In  this  section  the  results  of  the  analysis  are  presented.  The  analysis  procedure  is  first 
summarised,  followed  by  a  broad  look  at  some  general  aspects  related  to  successful  parameter 
estimation.  A  detailed  discussion  of  the  results  concludes  this  section. 

4.1  Analysis  Procedure 

As  indicated  in  previous  sections,  pitch  acceleration,  q,  is  required  for  corrections  to  normal 
accelerometer  readings  (Equation  (9))  and  angle  of  incidence,  «,  is  required  for  modelling  the 
non-linear  pitching  moment  and  Z-force  curves  (Equations  (3)  and  (4)).  Because  of  the  lack 
of  measured  data  for  these  two  quantities  a  two-stage  analysis  process  was  devised.  In  the  first 
stage  a  linear  model  was  used  and  q  corrections  to  normal  acceleration  measurements  were 
neglected.  The  resulting  match  against  measured  time  histories  of  pitch  rate  and  normal 
acceleration,  while  not  as  good  as  ultimately  obtained,  produced  calculated  values  for  q  and 
a.  These  were  then  used  to  make  the  required  correction  to  the  normal  acceleration  and  to 
set  up  the  full  non-linear  model  (Equation  (6)).  The  time  history  match  was  then  repeated  to 
complete  the  second  stage.  If  necessary,  a  further  iteration  could  be  made  based  on  the  calcu¬ 
lated  results  of  the  second  stage,  and  so  on.  However,  in  the  present  case  the  two-stage  process 
was  found  to  produce  a  sufficiently  good  match. 

4.2  General  Aspects 

In  Reference  10  a  number  of  conclusions  were  reached  relating  to  data  requirements  for 
successful  identification.  These  included  criteria  for  control  input  frequency,  data  sampling 
rates  and  record  length.  The  present  data  can  be  examined  in  the  light  of  those  criteria  and 
provide  some  experimental  validation  of  them.  Record  length,  T ,  in  all  cases  satisfied  the  require¬ 
ment  w„T  -■  14,  i.e.  about  two  cycles  of  oscillation  or  more.  The  criterion  on  sampling  rate 
suggested  in  Reference  10  was  w„Ai  <  014  and,  based  on  the  natural  frequencies,  <«„,  shown 
in  Figure  2  and  the  sampling  rate  of  40  per  second,  this  requirement  was  not  met  in  some  cases 
at  the  higher  values  of  o>n.  The  worst  cases  were,  as  is  to  be  expected,  at  the  higher  Mach  numbers 
and/or  lower  altitude.  The  results  of  Reference  10  further  indicated  that  fit  error  deteriorated 
fairly  gradually  as  ainA/  increased  so  that  some  relaxation  of  the  sampling  rate  criterion  can 
be  argued  provided  that  other  requirements  are  satisfied.  A  third  factor,  related  to  control  input 
application,  is  the  ratio  of  the  applied  doublet  frequency,  <or,  to  the  natural  frequency,  <«n- 
Reference  10  found  that  fit  error  deteriorated  rapidly  for  small  u>i/w„  and  below  a  certain  value, 
successful  identification  was  not  possible.  For  a  doublet  type  input  this  minimum  ratio  was 
approximately  0-4.  In  the  present  series  of  tests  the  pilot  had  difficulty  applying  a  sufficiently 
rapid  input  in  some  cases  due  to  the  high  control  forces  involved.  This  was  particularly  true  for 
the  low  altitude,  high  Mach  number  cases  where  problems  also  arose  with  the  sampling  rate 
criterion.  Examination  of  the  matched  results  showed  that  an  inadequate  match  could  be 
expected  for  wi/w„  below  0-4  to  0-45  even  when  w„A/  was  less  than  0  - 14.  On  the  other  hand  for 
larger  values  of  wt/wn  a  satisfactory  match  could  be  obtained  with  con  A /  as  high  as  018.  Because 
of  these  observations,  test  results  where  wt/u>n  was  less  than  0-4  have  been  rejected.  On  the 
other  hand  the  maximum  value  for<onA/  was  set  at  0- 18  provided  that  a>f/a>„  was  greater  than  0-45. 
Application  of  these  criteria  meant  that  all  results  at  10,000  ft  above  M  =  0-94  were  rejected 
and  some  high  Mach  number  results  at  30,000  ft.  were  also  rejected.  A  higher  sampling  rate  may 
have  saved  some  of  these  points  but  the  main  restriction  was  the  difficulty  encountered  in 
applying  a  sufficiently  rapid  control  input. 


Another  important  aspect  bearing  on  the  accuracy  of  the  results  is  the  maintenance  of 
constant  Mach  number  during  the  course  of  a  particular  test  case.  This  is  particularly  important 
in  the  transonic  regime  where  rapid  aerodynamic  changes  take  place.  Because  of  this,  a  careful 
examination  of  the  measured  time  history  of  indicated  airspeed  (IAS)  for  each  test  was  made 
and,  for  Mach  number  between  0-90  and  1  -05,  if  IAS  changed  by  more  than  1  %  in  the  course 
of  the  test,  the  results  were  rejected.  This  led  to  the  rejection  of  a  further  three  sets  of  results. 
An  exception  was  made  for  a  series  of  doublet  input  manoeuvres  at  M  approximately  equal  to 
0-95.  The  results  were  found  to  be  particularly  sensitive  to  speed  at  this  point.  A  slight  variation 
of  Mach  number  from  one  manoeuvre  to  the  next  produced  a  large  variation  in  the  results, 
thus  indicating  the  rapidity  with  which  some  derivatives,  in  particular  the  pitch  damping,  are 
changing.  This  will  be  discussed  further  in  Section  4.4. 


4.3  Classification  of  Results 

The  results  remaining  after  application  of  the  criteria  of  Section  4.2  are  summarised  in 
Table  3  which  presents,  in  columns,  the  true  Mach  number  (Mr)  together  with  the  non- 
dimensional  derivatives  and  “Cramer-Rao"  bounds  obtained  for  each  successfully  identified  case. 
Both  10,000  ft.  and  30,000  ft.  results  are  included.  The  identifying  numbers  at  the  top  of  each 
column  correspond  to  the  numbers  appearing  in  Figures  I  and  2.  Numbers  beginning  with  1  or  3 
refer  to  tests  at  30,000  ft  while  those  beginning  with  2  are  at  10,000  ft. 
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Two  examples  of  matched  results  obtained  are  shown  in  Figures  3  and  4.  Figure  3  corres¬ 
ponds  to  case  24  of  Table  3,  at  Afr  =  0-93  and  10,000  ft.  This  case  represents  a  marginally 
acceptable  match  with  the  ratio  wt /<o„  being  low  at  about  0-45  and  <o„A/  being  013.  The 
resolution  of  the  transducer  measuring  elevator  angle  was  not  as  good  as  may  be  generally 
expected  as  can  be  seen  in  Figure  3.  However,  reasonable  results  were  obtained  without  recourse 
to  further  processing.  Figure  4  shows  the  time  history  matches  for  case  310A  at  M  =  1-148  and 
30,000  ft.  Because  a  full  doublet  input  was  not  obtained  in  this  case  it  is  difficult  to  calculate  a 
value  for  However,  it  may  be  estimated  to  be  around  0-5  to  0-6,  a  quite  adequate  value. 
The  value  for  wn\t  is  somewhat  high  at  about  017  but,  nevertheless,  an  acceptable  match  was 
obtained  for  this  case. 

Figures  5  to  9  are  plots  of  the  results  summarised  in  Table  3.  In  addition,  the  a  priori  values 
of  each  of  the  non-dimensional  derivatives  have  been  plotted  in  these  figures  and,  for  reference, 
faired  estimates,  obtained  using  the  Digital  Datcom  Program11  are  also  shown  in  Figures  5,  7 
and  8. 


4.4  Detailed  Discussion 

4.4.1  Cramer-Rao  Bounds 

In  order  to  assess  the  level  of  confidence  to  be  placed  in  the  values  of  the  respective  derivatives 
extracted  it  is  helpful  to  examine  the  corresponding  “Cramer-Rao”  bounds.  Significant  reduction 
of  these  bounds  below  the  comparable  a  priori  weightings  indicates  the  presence  of  relevant 
information  in  the  measured  records  and  implies  a  reasonable  level  of  confidence  in  the 
extracted  value.  On  this  basis  it  is  clear  from  a  comparison  between  the  bounds  shown  in  Table  3 
and  the  a  priori  weightings  summarised  in  Table  2  that  a  fair  amount  of  confidence  can  be 
placed  in  the  pitching  moment  derivatives  (Figs.  5  to  7)  including  the  pitch  damping.  Of  the 
Z-force  derivatives,  only  the  C*a  results  (Fig.  8)  can  be  accepted  with  some  confidence  while 
the  results  for  CIte  indicate  very  little  information  for  this  derivative.  Consequently  its  extracted 
values  do  not  vary  significantly  from  the  a  priori  settings  (Fig.  9).  It  should  be  noted  that  intro¬ 
duction  of  a  larger  number  of  parameters  when  using  the  full  non-linear  system  of  Equation  (6) 
caused  some  increase  in  the  “Cramer-Rao”  bounds  as  compared  to  those  obtained  using  the 
simpler  linear  model.  This  was  particularly  true  of  the  pitch  stiffness  Cm*  which  is  most  affected 
by  the  non-linear  modelling.  “Cramer-Rao”  bounds  shown  in  Table  3  for  Cm*  are  up  to  twice 
as  large  as  those  produced  by  the  first  stage  linear  model,  but  are  still  considerably  smaller  than 
the  a  priori  weights. 

4.4.2  Pitch  Damping 

Figure  5  summarises  all  the  results  obtained  for  pitch  damping  at  both  10,000  ft  and 
30,000ft.  The  10,000ft  points  obtained  at  Afr  =  0-7  and  Afr  =  0-79  lie  slightly  below  the 
30,000  ft  data  indicating  a  possible  non-linearity  with  incidence.  However,  the  results  lie  within 
the  95%  confidence  interval  obtained  by  application  of  Student's  /-test  to  the  results  at  each 
Mach  number.  The  95%  confidence  intervals  have  been  evaluated  wherever  there  are  three  or 
more  data  points  at  approximately  the  same  Mach  number,  and  are  indicated  by  the  bars  in 
Figure  5.  In  general  the  scatter  of  the  results  can  be  regarded  as  good  with  the  relatively  large 
scatter  at  the  highest  Mach  number  attributable  to  the  fact  that  both  a>r/a>n  and  <unA/  are 
marginal  there.  The  results  at  M  —  0-95  are  exceptional  due  to  the  apparently  large  scatter. 

However,  these  results,  which  were  obtained  from  successive  doublet  inputs  at  a  single  test 
point  reflect  slight  variations  in  airspeed  from  one  input  to  the  next  and  indicate  the  very  rapid 
aerodynamic  changes  taking  place  at  this  Mach  number. 

In  general,  the  results  show  extracted  pitch  damping  values  to  be  quite  different  from  the 
a  priori  expectation.  At  subsonic  Mach  numbers  the  values  are  little  more  than  half  the  a  priori 
values,  and  rise  gradually  to  a  peak  at  about  Afr  =  0-94.  This  is  followed  by  a  rapid  drop  at 
around  Mi  =  0-95  and  a  recovery  to  reach  the  supersonic  value  at  about  Mi  =  0-96  to  0-97.  I 

The  supersonic  value  is  also  lower,  by  about  20%,  than  the  a  priori  value.  For  comparison,  the 
Digital  Datcom  results  appear  to  be  quite  close  to  the  extracted  values  at  subsonic  and  super¬ 
sonic  speeds  but  the  variation  at  transonic  Mach  numbers  is  quite  different. 
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4.4.3  Other  Derivatives 

Pitching  moment  due  to  elevator  deflection,  CmJe,  is  shown  in  Figure  6.  Although  the 
extracted  value  at  M  =  0  •  7  is  about  1 5%  lower  than  the  a  priori  value,  the  general  agreement 
at  subsonic  speeds  is  reasonably  good.  However,  for  Mach  numbers  above  1  the  flight  values 
are  consistently  about  20%  below  the  a  priori  (wind  tunnel)  values.  The  10,000  ft  result  appears 
to  be  somewhat  lower  than  the  30,000  ft  results  at  M  =  0-93  but  at  the  lower  Mach  numbers 
differences  between  10,000  ft  and  30,000  ft  do  not  appear  to  be  significant. 

Results  for  pitching  moment  due  to  incidence,  Cma ,  are  summarised  in  Figure  7.  At  subsonic 
Mach  numbers  the  low  flight  values  reflect,  to  some  extent,  the  non-linear  pitching  moment 
curves  (see  Reference  1  for  a  more  detailed  discussion).  The  linear  model,  at  stage  I  calculations, 
generally  produced  significantly  higher  values  of  Cma  at  subsonic  speeds  and  the  ratio  of  the  two 
can  be  seen  in  Figure  10.  Figure  10  also  shows  that  for  Mach  number  greater  than  about  1  both 
the  linear  and  non-linear  models  produce  substantially  the  same  value  for  Cma,  as  would  be 
expected  for  a  linear  variation  of  pitching  moment  with  incidence.  Returning  to  Figure  7,  the 
rise  of  Cma  at  transonic  M  appears  to  be  steeper  than  expected  from  the  wind  tunnel  results. 
This  is  consistent  with  the  similarly  rapid  variation  of  pitch  damping  (Fig.  5)  at  about 
Mt  =  0-95.  At  supersonic  speeds,  Figure  7  shows  the  extracted  flight  values  to  be  consistently 
about  20%  below  the  tunnel  results,  as  was  also  noted  in  Figure  6  for  Cmle.  The  Digital  Datcom 
curve  shown  in  Figure  7  is  higher  than  either  flight  or  wind  tunnel  values  throughout  the  Mach 
number  range. 

The  variation  of  Z-force  with  incidence  is  shown  in  Figure  8.  The  subsonic  values  of  C% 
measured  in  flight  appear  to  be  somewhat  lower  than  the  tunnel  values  as  was  noted  with  Cm* 
(Fig.  7).  However,  non-linearities  are  smaller  in  this  case  and  the  differences  in  C%  obtained 
between  the  linear  and  non-linear  models  are  generally  contained  within  the  scatter  band.  The 
transonic  flight  values  do  not  differ  significantly  from  the  tunnel  values  but,  as  with  the  pitching 
moment  derivatives,  the  supersonic  values  are  consistently  about  20%  below  those  obtained 
from  tunnel  data.  Digital  Datcom  estimates  lie  above  both  tunnel  and  flight  results  throughout 
the  Mach  number  range. 

The  fact  that  each  of  the  derivatives,  Cmte,  Cm*  and  Cr„  are,  at  supersonic  speeds,  substan¬ 
tially  less  than  the  respective  wind  tunnel  data  suggests  the  possibility  of  an  aerodynamic  cause 
common  to  all.  For  example,  the  differences  may  be  due  to  variations  in  boattail  geometry  between 
tunnel  model  and  full  scale  aircraft,  sting  interference  effects,  or  differences  in  flow  development 
particularly  about  the  elevators  and/or  wing-body  and  fin-body  junctions  due  to  Reynolds 
number  differences.  Further  detailed  investigation  would  be  required  to  resolve  these  questions. 

Finally,  the  results  for  Z-force  due  to  elevator  deflection,  Gt<,  are  summarised  in  Figure  9, 
which  shows  good  agreement  between  flight  and  a  priori  values.  However,  as  indicated  in  Section 
4.4.1,  the  flight  measurements  contained  very  little  information  on  this  derivative  and  little 
confidence  can  be  placed  in  the  results  extracted  from  flight.  The  a  priori  weighting  in  this  case 
prevents  the  results  from  straying  too  far  from  the  a  priori  values. 


5.  CONCLUDING  REMARKS 

A  flight  investigation  has  been  undertaken  to  determine  the  pitch  damping  of  a  60-degree 
delta  wing  aircraft  in  the  Mach  number  range  from  M  =  0-7  to  M  =  1-2.  The  flight  tests,  includ¬ 
ing  instrumentation  and  data  recording,  were  conducted  by  ARDU  and  the  data  provided  to 
ARL  for  analysis.  Apart  from  pitch  damping,  results  were  also  obtained  for  pitching  moment 
and  Z-force  derivatives  with  respect  to  incidence  and  elevator  deflection.  In  addition,  further 
information  relating  to  the  effects  of  sampling  rate  and  control  input  on  the  identification 
process  was  obtained.  This  will  have  a  bearing  on  future  flight  test  planning. 

A  modified  Newton-Raphson  parameter  estimation  procedure  was  used  to  analyse  flight 
records  of  normal  acceleration  and  pitch  rate  obtained  in  response  to  an  elevator  input.  The 
mathematical  model  describing  the  aircraft  motion  made  provision  for  non-linearities  in  pitching 
moment  and  Z-force  with  incidence,  and  also  took  account  of  a  time  lag  in  the  pitch  rate 
measurements. 

The  results  for  pitch  damping  were  significantly  different  throughout  the  Mach  number 
range  from  the  values  previously  assumed.  The  new  values  are  up  to  50%  less  at  subsonic  speeds 
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and  about  20%  less  at  supersonic  speeds  with  a  rapid  transition  taking  place  at  about  M  -  0  95. 
Of  the  other  derivatives  extracted,  pitching  moment  changes  due  to  incidence  and  elevator 
deflection  and  Z-force  change  due  to  incidence  all  had  values  about  20%  below  the  respective 
wind  tunnel  results  at  Mach  numbers  above  about  I,  while  showing  closer  agreement  at  lower 
Mach  numbers.  Further  investigation  would  be  required  to  isolate  the  causes  of  the  differences 
at  supersonic  speeds.  The  final  derivative,  for  Z-force  change  with  elevator  deflection,  could  not 
be  extracted  with  confidence  due  to  insufficient  information  contained  in  the  measurements. 

In  summary,  the  parameter  estimation  procedure  has  proved  to  be  an  effective  analysis 
method  even  when  only  a  modest  amount  of  flight  data,  obtained  using  a  relatively  simple 
instrumentation  package,  has  been  available.  Most  derivatives  have  been  successfully  extracted 
with  relatively  small  amounts  of  scatter  and  trends  have  been  clearly  identified. 
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NOTATION 


Matrix  of  stability  parameters.  Equation  ( I ) 

Measurement  bias  vector 

Matrix  of  control  parameters,  Equation  (I) 

Unknown  parameter  vector 
Mean  aerodynamic  chord 

Non-dimensional  pitching  moment  coefficient,  M  lpye2&c 
Non-dimensional  Z-force  coefficient,  Z/JptVS 
Non-dimensional  moment  derivative  w.r.t.  q ,  SCm/J(qf/2k'e) 
Non-dimensional  moment  derivative  w.r.t.  a,  JCm/Ja 
Non-dimensional  moment  derivative  w.r.t.  a,  iCm/Z(af/2Ve) 
Non-dimensional  moment  derivative  w.r.t.  &e,  iCm;iSe 
Non-dimensional  Z-force  derivative  w.r.t.  a,  ZCt/ba 
Non-dimensional  Z-force  derivative  w.r.t.  6,,  JCj/iS, 

Matrix  of  stability  parameters.  Equation  (2) 

Gravitational  acceleration 

Matrix  of  control  parameters.  Equation  (2),  or  gain.  Equation  (S) 
Moment  of  inertia  in  pitch,  kg-m2 
Mass,  kg 

Mach  number  or  pitching  moment 
Equation  error  bias  term.  Equation  (6) 

Pitching  moment  derivative  w.r.t.  q 
Pitching  moment  derivative  w.r.t.  * 

Pitching  moment  derivative  w.r.t  a2 
Pitching  moment  derivative  w.r.t  «a 
Pitching  moment  derivative  w.r.t  a 
Pitching  moment  derivative  w.r.t.  8, 

Normal  acceleration  in  g  units 
Pitch  rate,  rad/s 
Reference  wing  area,  m2 
Time,  sec 

Length  of  measured  responses,  sec 
Vector  of  control  input  variables 
Resultant  airspeed,  knots 


x  Vector  of  state  variables 

x  Body  axis  co-ordinate  in  forward  direction 

y  Calculated  response  vector,  Equation  (2) 

z  Body  axis  co-ordinate  in  down  direction 

Z  Force  in  z-direction,  N 

Zo  Equation  error  bias  term,  Equation  (6) 

Z»  Z-force  derivative  w.r.t.  q 

Z„  Z-force  derivative  w.r.t.  * 

Z„ i  Z-force  derivative  w.r.t. 

Z„»  Z-force  derivative  w.r.t.  a* 

Z,r  Z- force  derivative  w.r.t.  S« 

a  Angle  of  incidence  increment,  rad 

A  Increment 

8e  Elevator  angle  increment,  positive  t.e.  down,  rad 

{  Damping  ratio 

0  Pitch  attitude  increment 

p  Air  density,  kg/ms 

a  Standard  deviation 

r  Time  constant,  sec 

tut  Doublet  forcing  frequency,  rad/s 

a>m  Natural  undamped  frequency,  rad/s 


Subscripts 

a  Normal  accelerometer 

e  Trim  or  equilibrium  state 

m  Measured  value 

T  True  value 
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FIG.  1 :  INCIDENCE  RANGE  OF  TESTS— SUMMARY 


2.0 


Yw 


n  » 


Hz 


1.8 


1.6 


x 


X 


1.4 


©  Run  1 1...  ...  , 

A  Run  3  ]  30,000  ft. 

□  Run  2— 10,000  ft 

x  Run  2— 10,000  ft 
(Based  on  A  Priori  estimates) 


x 


1.2 


1.0 


0.8 


A  312  j 

&311 

©111 

39^38 


X 


818 
37 

©017 
1 7 


A  36 


0  23 
024 


0.6 


0.4 


0.2 


15©  E 
16%  3! 

4 


022 

35 

34 

14 
33  Q21 


<$32 

”<*>31 


I _ _  i 

-  i»... 

0  0.2  0.4  0.6  0.8  1.0 


FIG.  2:  FREQUENCY  AND  DAMPING  OF  SHORT  PERIOD-SUMMARY 
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FIG.  4:  MATCHED  RESULTS— CASE  31 OA,  M  =  1-148,  30,000  ft. 
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FIG.  5:  VARIATION  OF  PITCH  DAMPING  WITH  MACH  NUMBER 


m5e  Per  degree 


0  Run  1 
A  Run  3 1 
□  Run  2— 10,000  ft. 


30,000  ft. 


0$ 


4D 

str¬ 


op 


X- 


„-x 
- - 


40 


\ 


■  A  Priori 


a 


% 


'A 

A 


G& 


\ 


-V 


0.7 


0.8  0.9  1.0  1.1 


1.2  1.3 


FIG.  6:  PITCHING  MOMENT  DUE  TO  ELEVATOR  vs  MACH  NUMBER 
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FIG.  8:  Z-FORCE  DUE  TO  INCIDENCE  vs  MACH  NUMBER 
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FIG.  10:  RATIO  OF  LINEAR  AND  NON-LINEAR  MODEL  RESULTS  FOR  C, 
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